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_ Analytic Solutions to Several Optimum 
Orbit Transfer Problems 


Ht. Munick, R. McGill and G. E. Taylor 


bstract 


The analytic solution is given for the absolute minimum 
laracteristic-velocity path from a point on an elliptical 
bit to a nonintersecting circular orbit. The orbits are con- 
dered coplanar and two impulses are used. Existence proofs 
re given followed by proofs that the absolute minimum 
tisfies both necessary and sufficient conditions. 

More generally the complete analytic solution to the 
riminal-to-terminal problem as first formulated by Vargo 
fis given. 

An interesting consequence of the analytic solution is the 
roof that in optimal transfer the impulses are not in general 
yplied tangentially. Impulses are applied tangentially only 
; terminals having zero radial velocity. 

Each point on an elliptical orbit corresponds to an arrival 
unt of an optimum path from a lower energy circular orbit. 
he Hohmann path arriving at apogee is rigorously shown to 
> the optimum of the infinity of optimal transfer paths, for a 
ide class of orbits. 


ymbols 


= normal component of velocity 
= radial component of velocity 
= nondimensional normal component of velocity 
just after first impulse 
= nondimensional radial component of velocity 
just after first impulse 
= initial radial distance 
= final radial distance 
= distance ratio = a/8 
= nondimensional characteristic velocity 
‘= characteristic velocity of transfer 
= gravitational constant 
= nondimensional normal component of arrival 
velocity 
= nondimensional radial component of arrival 
velocity 
bscripts as defined in text 


*Presented at the 11th International Astronautical 
meress, Stockholm, Sweden, August 14-20, 1960. 

+ Grumman Aircraft Engineering Corporation, Bethpage, 
sw York. 

1 Numbers in brackets indicate references at end of paper. 


Introduction 


The problem solved in this paper is that of finding 
the analytic expressions for the minimum characteris- 
tic-velocity path (hereafter called the optimum path) 
for a class of orbital transfers. The optimum path solu- 
tions are derived for transfer between two terminals, 
and between a terminal and a ‘‘point.’’ For the purposes 
of this paper we define a terminal as a locus specified 
by a radial distance and a velocity vector, i.e., a posi- 
tion on an ellipse with the line of apsides unspecified. 
A “point” is defined by a radius vector and a velocity 
vector. A “point”? then, is a specific position on an 
ellipse whose line of apsides zs specified. The class of 
orbits considered are those for which the major apsis 
of one terminal is less than the minor apsis of the other 
terminal. This is tantamount to saying that a minimum 
of two impulses is required to accomplish the transfer. 

The formulation given herein is that of Ref. [2] 
which was pioneered by Lawden. 

The assumptions of the paper are: 

1) Inverse square force field 

2) Two-body equations 

3) All orbits are coplanar 

4) Impulsive thrusting 

5) Two impulses are used. 


Definition of Problem 

Consider the problem of optimum transfer of a 
space vehicle between two terminals. Let (wo, vo) be 
the velocity components of the first terminal at a dis- 
tance a from the attractive center. Let (wr, vr) be the 
velocity components of the second terminal at a dis- 
tance 6 from the attractive center. At the first terminal 
an impulse is applied resulting in new velocity compo- 
nents (%4, v1). The space vehicle then goes into a 
transfer orbit, arriving at the second terminal with 
velocity components (wz, v2). Upon arrival at the 
second terminal a second impulse is applied adjusting 
the arrival velocity to the desired terminal velocity 
(uy, Ur). For such a maneuver the characteristic- 
velocity is given by 


AX = a/ (uu — UW)? + (V1 = Uo)” 
ap VJ/ (us — ur)? + (ve — vr)*. (1) 


From conservation of angular momentum we have 


U1 » (2) 


UU = 


BIR 
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and from the conservation of total energy 


ih 1 
oo =[1— S] ut + ott [5 - =p (3) 


Introduce dimensionless parameters x, y by: 


maa 4/*, naug/ts (4) 
a a 


and in all other cases: 
D, = Yi / ne (5) 
a 


w= g/t 
a 


Let the distance ratio be denoted by 


=<, (6) 
"ee 
Using (2), (4), (5), and (6), the characteristic- 
velocity becomes 
Naess we [VW (a = a)? + (y — Yo)? (7) 


+ V(rx — ar)? + (ye — yr)’. 
And Eq. (3) becomes 


Se iy 2a). (8) 
From (7) and (8) we get 
ot | 4/8 = Vea EU we (9) 


+ V(rx — tr)? + (y2 — yr). 


The sign of the radical 4/(1 — r?)a®. + y? + 2(r — 1) 
as obtained from (8), must be taken the same as y; . 
This is dictated by the fact that we are minimizing 
characteristic-velocity. 

Let A(x, y) denote the nondimensional characteristic- 
velocity as given by the right-hand side of (9). From 
(8) we observe the constraint relationship that 


(1 — P)a? + y + 2(r — 1) ZO. (10) 


If (10) is violated an imaginary radial velocity results. 
The physical interpretation of this is that the vehicle 
cannot achieve a distance 6 from the focus; hence a 
transfer in this case is impossible. 

From (8) and (9), we can now define the mathe- 
matical problem of interest to be that of minimizing 


Ma, y) =V (x — ao)? + (y¥ — yo)? 
oh Vitae tr)? + (yy — Yr)? (11) 


subject to 


(l—P)a? + y? + 2(r — Dee: (2) 


with 
CO) <p < il. (GIES) 
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Solution to Problem 


In [2] it is rigorously proved that the functac 
A(x, y) given by (11), will assume a relative minimt 
at an interior point bounded by the closed cury 

= Oand y = yr. Since A(z, y) is a differentiah 
eke in such a region, a necessary condition for th 
solution (x, y) of the minimum problem is given by: 


It is also shown in [2] that using the two necess 
conditions given by (14) we establish for y. # 0 


(1 = Par? + 2(r = DIK? : 
— [2ye{(1 — r°)rara + 2(r — 1)})K ‘ 
+ [0 — P)ryee? + 2 — Dyel=0 @ 
where 
Rea Pie ae oe 
Introducing the transformation yo’ = (y2/yr) m 
(15), and ies tee we get 
(1 + re? — Qly.” — 2y [1 + r)rert — 2] Rt 


+ [1 + r)r2? — 2] = 0. 
Re-arranging terms 


(1 + r)r]e” — Qa[(1 + r) rary’ 


; (Vv 
+ [1 + r)crys — 2(yo' — 11 = 
Dividing through by (1 + r) we get 
: 2 
re — Qrapyo'e + E yo + (yo’ —1) = an E 
ied ear (1 
XL SY = i Bs as —— = 0 
x E Yo (yo 1) Vj ae 
Factoring (18) yields 
[r= — {aru + (ue! — 1) 4/72} 
To Sm yr! Yo" = 
a 1l+r al 


x [re — fae 5 ee Wor ahl-9 


From (19) we conclude that the minimum f 
(av, y) must satisfy one of two possibilities: 


tr Yo’ + (ye! — 1) Vem (2 
ir Yo — (ye! — 1) oe 


Combining terms in (20), and using yo’ = 


> 
cs 


or 


= 
8 
| 


(Y>/ Yr) 


be ee 2) (2. 2 
Pen 1 (ar + 4/72). (2 


Further simplification of (22) yields 


» 
ss 
UP Se he Yr 

TL — Lp 9 (23) 
4 _—_———_ 
“4 Wa = / any, 
sing (23) and the expressions for the derivatives 
h/dx), (A/dy) given in [2] we conclude 
* 


a Y— Yo _ Yo 
w — Xo THT Rs (24) 
mtr 4/2 


Hq. (24) gives y as a linear function of a, and setting 
us linear function of x into (23) yields a quadratic 
juation in x whose coefficients depend on the given 
ita. We then solve analytically and determine two 
cact values of x. 

Using the second possibility given by (21) and pro- 
eding as above we obtain 


Ure E Yr 
TL — Xp / 9 (25) 
tr — — 
1l+pr 
¥Y — Yo Yo 


Similarly the pair of Eqs. (25), (26) can be solved 
ialytically yielding two more exact values of x. Using 
44), (26) the corresponding y values are determined. 
We have now determined analytically four solutions 
*, y) to the minimum problem. Using these four 
ilues of (x, y) there correspond four values of \. The 
nallest of these four values of ) is the absolute minomum 
r \. This then gives the complete solution to the 
rminal-to-terminal problem first formulated by 
argo [1]. It will be shown subsequently that under 
rtain orbital conditions three of the four values of \ 
n be rejected analytically. We next study geometri- 
lly the optimum transfer path given by the above 
ialytical expressions. 

Replacing rz by x2 in (23), we have the pair of equa- 
m from (23), (24) 


Y2 — UF - UF 
to — Up Lee 
1 y 3 
ae (27) 
Y— Yo Yo 


Rqs. (27) tells us that only in the case where the 
st or second terminals have radial velocities ap- 
oaching zero will the transfer orbit leave or arrive 
ngentially. In all other cases the transfer orbit does 
t leave the first terminal tangentially nor arrive at 
e second terminal tangentially. In the interesting 
se of transfer between a circular orbit and a point on 
nonintersecting elliptic orbit, the transfer path is 
igential to the circular orbit, but is not tangent to the 
iptical orbit (see Figure 1). 


Nontangential Arrival 


Tangential Departure 


Fig. 1. Optimum Circle-to-Ellipse Transfer 


Eqs. (25), (26) give the only remaining possibility, 
and the same tangency conditions hold for this case. 
We next study the case of transfer between an elliptic 
and a circular orbit. 


Ellipse—Circle Transfer 


For the case where the second terminal corresponds 
to a circular orbit we have 


ary eae ea a 
B a B a 


vy = 


Nondimensionalizing we get 
1/2 
tp =r, Yr = 0. 
Since (15) does not hold for y2 = 0 we first consider 
transfer to an almost circular orbit 
ee nae Un =e 
where e¢ is an arbitrarily small radial velocity. In this 
case the annular region wherein the solution lies is 
bounded by two ellipses arbitrarily close to one another. 
Also the second possibility given by (25), (26) can be 
eliminated for positive yr . It is shown in [2] that y2 is 
less than yr , and rz is less than wp . That is, the trans- 
fer orbit is itself an ellipse. Therefore the left-hand side 
of (25) is strictly positive. Since 7 must le between 
zero and one we conclude 


ae 1/2 We ZO 
tr Se =— . 
ME l+pr ih + r 


Therefore (25) is physically unrealizable. Eq. (23) 
for arbitrarily small radial velocity yr becomes 


Ya = € € 


eo le es 4/ Ors (28) 
; IL Se % 
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Since € can be chosen arbitrarily small we conclude 
from (28) that in transfer to a circular orbit 2 becomes 
arbitrarily small. We then say 


i = ate eer a an Oe) 
and from (24) 


eS ante, (30) 
xv — Xo 
where 
tan go = Yo 


yee (31) 
One ag 


To determine the possible roots we first introduce 


the following 
r /; 2 , tan go. (32) 


Let A = tan ¢o 
B 
Using these values of A, B we get from (29) 


I 


Yo — Lo tan Co. = 


(1 — 7 + A*)x’ + 2ABr + [B’ — 211 — r)] = 0. (33) 


Solving for x 


ants —AB 
Tae A 
22 2 2 (34) 
oe Ee ete EEE Sele ON 
1-—r4+ A? 


From (32), (34) and some simplification we get 


2 y —_ 
Daas maa! r) V2(1 +17) sec go (35) 


ge 


— sec? go 


Since the denominator in (35) is strictly negative we 
choose the negative possibility for the numerator to 
insure that x is given as a positive quantity. That x 
must be positive is determined by the fact that in this 
case v and wp have the same sign (both orbital rota- 
tions taken in same sense). lor the retrograde case 29 
and x» would have opposite signs and x could be nega- 
tive. This gives 


& ater 5 meee, ve pee ( 


ee sec? po, 


36) 


In general the nondimensional characteristic-velocity 
can be written as 


= n= (em) /1 +(' = 
(’ = | 
Eire <= oe) 1 ee (eaey. 
V2 — Up 


For the ellipse to circle case using (23), (24), (36) and 
(37) the absolute minimum for ) is given by 


(37) 
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\absolute 
minimum 


2 
1/2 
é +o 
— (mtr os) +0 


For the special case of a = 0, i.e., the initial veloc 
all radial, this result reduces to that given by Lawd 
[3]. Also for the case wherein yo = 0, circle-to-cire 
Eq. (38) reduces to the Hohmann result which 
hereby clearly shown to be an absolute minimum. 

An analogous result may be obtained for the case 
circle-to-ellipse transfer. 


Optimum Apogee Arrival Paths 


ee ee 


In [2] it is demonstrated that for optimum circle-t 
ellipse transfer the first impulse must be appli 
tangentially. For this case the nondimensional chara 
teristic-velocity is found to be 


A=a—14+WA(re — ar)? + (Yo — Yr): (3 
From (39) 
N2a—1+4r —172, (4 


the nondimensional velocity xr , in terms of the orbs 
elements, is given by 


np sey / ee (4 
lp 


where: er = eccentricity of elliptic orbit 


I 


rp = distance ratio at perigee. 


Since departure from the circular orbit is tangenti 
(y = 0) it follows that 


2 
aes ee 4 
reas 


This may be deduced from a study of the annul; 
region containing the solution (see Ref. [2]). Combinin 
(39), (40) and (41): 


ie hey l+e 
IN es aS —r ; 4 : 
= ery rr 4/ 1. (4 


The least value of the derivative with respect to r 
the right-hand side of (43) is 


If the least value of this derivative is positive then ; 
values of the derivative are positive. Therefore t 
right-hand side of (43) is monotonically increasing ir 
provided 


a 
45° 


se 
IIA 


(4 


is shown in [2] that at apogee or perigee the equality 
(43) holds. For the class of orbits given by (44) the 
if ht-hand side of (43) assumes its minimum value at 
he smallest r, i.e., at apogee. Since from (48) it is 
mpossible for \ to assume a smaller value than the 
ainimum of the right-hand side it follows that \ as- 
umes its least value at the smallest value of r, which 
ccurs at apogee. 
“Starting from the circular orbit, the methods given 
i this paper give the optimum path to any given point 
n the elliptical orbit. Therefore to every point on the 
liptical orbit there corresponds an optimum transfer 
ath. Of the infinity of optimal paths there is an 
ptimum transfer path. It is rigorously demonstrated 
bove that for the class of orbits defined by (44) the 


Hohmann path arriving at apogee is such a least fuel 
path. 
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Orbital Patterns for Satellite Systems 


Louis G. Vargo* 


Abstract 


Groupings of satellites are considered as a means of per- 
forming the equipment-carrying function in several types of 
advanced systems. These groupings and their associated 
orbital arrangements or patterns are analyzed on the basis 
of their ability to continuously cover a large part of the 
Earth’s surface. Optimal patterns are determined which 
minimize the total number of satellites required. Some 
operational features of these patterns are noted with the 
result that no fundamental technological barriers appear to 
exist in their realization. It is also concluded that determinate 
patterns as opposed to those involving probabilistic coverage 
are preferable. 


Introduction 


Several types of proposed space systems demand a 
continuous or nearly-continuous capability to gather 
information over large areas of the Earth’s surface and, 
in some cases, transfer this information to other widely 
separated locations. Communication systems using 
line-of-sight transmission are typical. These systems 
have been investigated by Haviland! and Vilbig? with 
some treatment of the orbital problems encountered 
when terrestrial satellites are employed as equipment 
carriers. The central purpose of the present paper is to 
examine the orbital aspects of a broad class of systems 
which have the above-mentioned coverage require- 
ments. This class may also include navigation, metero- 
logical, and military reconnaissance, surveillance and 
bombing objectives. 


Instantaneous Coverage Areas 


In general, the geographic area from which a satellite 
can obtain data (or, conversely, the locations which can 
receive information from the satellite) at any instant of 
time is a function of the position and velocity of the 
satellite relative to the Earth. For most cases, this area 
is quite insensitive to satellite speed and direction be- 
cause of the great relative speed of the band of electro- 
magnetic radiation used for sensing and transmission. 
Systems employing some type of scan device require a 
closer determination of the velocity-area dependence. 
Scan rates of the order of one cycle/sec. or greater repre- 
sent current practice. This value and the greatest angu- 
lar speed associated with a low-altitude circular orbit 
imply a displacement of the sub-satellite point of about 
five nautical miles per cycle. We will show that this 
displacement is at least two orders of magnitude down 
from the linear dimensions of a typical area seen by the 

* Research Engineer, Astrodynamics Department, Aeronu- 


tronic, A division of Ford Motor Company, Newport Beach, 
California. 
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equipment during any one cycle. On this basis one ma 
ignore the velocity effect in the preliminary desig 
phase of the orbital system and associate an area of th 
Earth’s surface with only the instantaneous position ¢ 
the satellite. This area will be referred to as the instar 
taneous coverage area (ICA). It contains all geographi 
locations at which the system purpose is accomplishe 
by a single satellite at one instant. ; 

Many ICA’s exhibit radial symmetry about the suk 
satellite point. This may result from the use of fre 
rotational scanning for purposes of injecting a constan 
angular momentum vector into the attitude contre 
dynamics of the vehicle. The dynamical analysis is thu 
simplified when compared with that required for sca 
techniques involving reciprocative motion. For a give 
scan rate, peak component accelerations and problem 
related to bearing design are alleviated. The rotatione 
axis may then be directed along the local vertical anc 
if slant range or incidence angle limitations are imposec 
these conditions can be met uniformly in all azimuthe 
directions. Fig. 1 shows the geometric situation alon 
with definitions of the pertinent quantities. 


— 


bi! 
tH 
a9) ef 


ITH 


“ie 


SATELLITE 


Fia. 1. Circular instantaneous coverage areas 
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Fig. 2, Correspondence between circular ICA size and satel- 
ite distance for various incidence angle and slant range limits. 


The spherical surface of the Earth implies circular 
[CA’s for the above type of satellite equipment. This 
‘orm of ICA is shared by systems employing fixed an- 
fennas which are essentially omnidirectional. Line-of- 
sight range and incidence angle restrictions arise on 
ircles about the sub-satellite points. If these restric- 
ions are not present, the ICA is bounded by the local 
yorizon and, again, is circular. The curves of Fig. 2 
we useful in determining trade-offs between equipment 
yerformance factors, satellite altitude and the ICA 
ize. We plot the half-central angle ¢ vs. radial distance 
-for various maximum incidence angles with maximum 
Jant range curves superimposed. For a given type of 
quipment, the orbital patterns for some overall cover- 
we requirement evolve by consideraion of the ¢ — r 
‘combinations within the slant range and incidence 
gle limits.* The criteria and procedures for pattern 
election from these ICA-altitude combinations are 
liseussed in subsequent sections of the paper. 


* Other limits peculiar to the system may arise. The emitted 
R radiation of ascending rocket vehicles may be the informa- 
ion sought in a military surveillance system.* Here a “ground 
lutter’”’ or background noise problem exists, so that ICA’s take 
n annular form. The equipment is blind within the inner circle. 
“his circle must not be so large that a vehicle may escape detec- 
ion, and consequently the outer circle is delimited. 


One may imagine other than circular ICA’s produced 
by various equipment designs and operational tech- 
niques. It is felt they are of secondary importance be- 
cause of the many factors leading to the circular ICA. 
However, if noncircular cases do arise, the following 
discussion may be taken as a model for their inclusion. 


Pattern Criteria and Synthesis 


We may now address the problem of satisfying a sys- 
tem coverage requirement with an assemblage of circu- 
lar ICA’s. Each ICA shows a time dependence governed 
by the orbital motion of its associated satellite. For an 
example, assume that the requirement is continuous in 
time and global in extent. That is, every location on the 
Earth must lie within an ICA at every instant after 
completion of the pattern. As a starting point, we as- 
sume the orbital motion to be purely Keplerian and re- 
gard all non-inverse-square forces as possible disturbing 
causes. If the required coverage area is symmetric with 
respect to the Earth’s center, then ICA’s of constant 
diameter for all satellites and within the orbit of each 
satellite should be employed. This is implied by the fact 
that the locus of sub-satellite points does not, in general, 
retrace itself. (The exceptional cases given in a recent 
paper‘ have a low probability of occurence.) Thus, orbits 
have no preferred locations within the latitudinal limits 
set by their inclinations over a time interval which is 
long compared to the orbital period. Since the area be- 
tween these limits is also centrally-symmetric, a “prin- 
ciple of uniformity”? may be invoked. This leads to the 
consideration of patterns made up of circular orbits of 
equal altitude. It should be emphasized that this ap- 
plies only to the case of a centrally-symmetric coverage 
requirement. 

A global coverage requirement may be satisfied by 
many patterns. Selection from these patterns involves 
consideration of diverse criteria. An attempt will be 
made to state the basic criteria and determine which 
should be used in the initial pattern selection and which 
should be assigned to more advanced phases of the sys- 
tem design process. 

(a) Attainment Effort (costs) 
number of satellites in the pattern, the equip- 
ment type and weight in each satellite, and the 
orbital altitude. 

(b) Maintenance Effort (costs)—dependent on many 
of the factors in the attainment effort but with a 
strong statistical basis. 

(c) Stability—essentially measured by the time in- 
terval from an initially satisfactory pattern to 
violation of the coverage requirements due to 
pattern perturbations—not associated with 
equipment reliability. 

(d) Compatibility with present and future launching 
facilities and tracking ranges. 

(e) Flexibility—the relative — effort 
adapting to other system requirements. 


a function of the 


involved in 
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One quickly concludes from these criteria that (a) is 
the only one suitable for pattern initial selection, Items 
(b) and (c), though essentially quantitative in nature, 
require an assumed set of orbits to obtain numerical 
values; items (d) and (e) are essentially qualitative. It 
will be shown that use of (a) allows for a direct synthesis 
of optimal patterns so that design charts may be gen- 
erated. These then serve as a basis for the consideration 
of the other criteria. 

Pattern selection from a fixed ICA diameter implies 
prior specification of equipment coverage characteristics 
and operating altitude. The remaining attainment 
effort parameter, viz. number of satellites in the pattern, 
may then be minimized. The costs related to these mini- 
mum numbers can be compared with the costs involved 
in developing equipment with various ICA capability. 

A “ring” is defined to be a group of equally-spaced 
satellites whose orbits differ only in the time of nodal 
passage. Fig. 3 shows the relative positions of the ICA’s 
for such a group composed of m satellites. The ICA 
centers are separated by a central angle of 2 /m 
radians. If m is large enough so that this spacing is less 
than the ICA diameter 2 ¢, a continuous strip of width 
2w will be covered. The particular rig shown makes an 
angle a with the equator where a is also the orbital in- 
clination of all satellites in the ring. Q is the nodal 
longitude measured from some reference meridian. 
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Fra. 3. Ring formation from equally-spaced ICA’s 
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Fra. 4. Geometry of circular ICA overlap 


The ring half-width is a function only of g and m. Br 
cause of the curvature of the Earth’s surface, two ci 
cular ICA boundaries, projected on a plane normal 1 
the local vertical at their line-of-intersection (A-A 1 
lig. 4), will be ellipses with equations: 


2 2 
Es y 

eee sar tS 1 

FR? sin? ¢ cos? 2r/m R? sin? ¢ . ( 
(2 — 2z)? y 

+ lst =a) (2 


R? sin? ¢ cos? r/m sR? sin’ ¢ 


R is the radius of the Earth. From the figure, d = 
R cos ¢; therefore, = dsinr/m = R cos¢ sinz/m. Th 
ellipses intersect at the points, «7 = = R cos¢ si 
r/m, Yy = +9 = +R (sin’e—cos’g tan? x/m)', so the 


w = sin (sing —cos’¢ tan? 2/m)? 
= sin! 8 fe,m). @ 


Now rings may be assembled to effect the require 
global coverage in two ways. First, several rings of th 
same inclination but different nodal longitudes may I 
combined to extend coverage longitudinally. Secon¢ 
the nodal longitudes may be held fixed, but their i 
clinations varied to provide greater latitudinal coverag: 
A pattern which obtains global coverage by the fir 
method we class as an L-type. If the second is employet 
a d-type pattern results. It will be shown that rath 
small differences exist between the minimal number | 
satellites for each type; certain operational consider 
tions dominate any selection. 


To determine the optimal patterns for each type, 
he critical region of least overcoverage must be 
estigated. In the case of the L-type, a ring inclina- 
n of a = r/2 — wis prescribed. The number n, of 
rings required depends on the sector of the 
quator corresponding to each ring. From Fig. 5, 


= sin w 
Site: = tan w. (4) 
sin @ 


Bing N, be the minimum total number of satellites 
equired and noting that 2 na = for complete coverage 
round the equator, 


L< 


wd 


> Mm 7 
P. a= Te : : ; 
a 2 sin tan sin“ 6(y, m) 


a 


(5) 


Nearly this lower bound is defined only at values of 
and mfor which tan sin-'8 < 1, ie. w S$ 45°. When 
p= 45°, eq. (4) implies « = 90°, and the two types of 
atterns coalesce. 

~The critical region for A-type patterns is along the 
elestial meridian of anti-nodes common to all rings. 
gain, if n is the number of rings, 2nw = 7 gives 
lobal coverage, and 


m 1 
Ny = 2=—— : j 
Be ES gin B(y, m) (6) 


Fra. 5. Local geometry of equator for L-type patterns 
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Fia. 6. L-type pattern optimization 


Comparing eqs. (5) and (6), we see that the lower 
bound for N, is less than the bound for Ny. The full 
mathematical formulation of the problem of determin- 
ing optimal patterns must include the integer property 
of m, n, N, and Ny. Relations between m and n must 
also be satisfied, viz. 


n= fi (7) 


2 sin tan sin 6’ 


IIV 


and 


a (8) 


2a Bo 


IV 


for L- and )-types respectively. A straight-forward 
solution process consists of exploring values of NV; and 
Ny for integer combinations of m and n which satisfy 
eqs. (5-8) regarding ¢ as a parameter. When this is 
done the dashed curves of Figs. 6 and 7 result. It is 
interesting to compare these with curves obtained by 
assuming the variables to be continuous. These are 
shown as solid lines in the figures. The lowest point of 
the dashed curves gives the optimal pattern variables. 
In some instances these variables are not unique. 
When ¢ = 40° for example, an N, of 18 results from 
mn = 6 X 3 and 9 X 2. Both arrangements meet the 
overall coverage requirement with the same minimum 
total number of satellites. Fig. 8 gives N; and Ny as 
they depend on ¢ in the range 10° <¢ 40°. Because 
of the relatively large magnitudes of N;, and N) in this 
range, a smooth curve through several computed points 
allows intermediate values to be estimated quite 
accurately. When ¢ is greater than about 40°, however, 
the step-wise character of the minimal numbers needs 
exact representation. 
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Fia. 7. A-type pattern optimization 


The process of integerizing the solutions produces 
some ring overlap at the critical regions even for optimal 
values of the variables. This degree of overlap depends 
on the parameter ¢ as would be expected. Subsequent 
sections of the paper will consider the manner in which 
this overcoverage effects pattern stability, attainment 
accuracy, ete. 


TOTAL NUMBER OF SATELLITES IN OPTIMAL PATTERNS 


Np AND Ny, 
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Fra. 8. Minimal number of satellites vs. ICA size 
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TABLE 1 


Optimal Pattern Variables for Continuous Global Covera 


g Circular w, Ring m., Number Neely N,T 
A Half- Patt ipa * n, Number) numb 
central Angle, Type ‘ BePae es Rae of Rings, | “Satellit 
egs. 
10 iv y 15 13 
7 24 14 
: 8 28 12 
20 ic 15 13 6 
14 12 i; 
"e 16 14 <6 
30 iy 22 9 4 36 
N 24 10 4 40 § 
40 28 6 3 138 
¥ fe 9 2 } 
d 32 7 3 21 
50 i 37 5 2 10 
d 37 5 3 15 
60 L 44 4 2 Sa 
it 52 5 2 10. 
6 
70 L\ ie 3 2 
: 
80 iS f 3 2 6 
xi 69 
2 
90 a a y 1 


Table 1 gives the values of the pattern variables f 
minimum total number of satellites for g¢ = 10°, 2( 
30°,---, 90°. Between ¢ = 60° and 70°, w for optim 
L-type becomes greater than 45°. Although the ts 
pattern classes lose their identity at this point, t 
manner in which the intrinsic overcoverage is distri 
uted retains their features. When, for example, ¢ 
80°, we may place both rings at an inclination a = 9 
— w = 21° to just cover the poles and obtain douk 
coverage in the equatorial region. On the other han 
the inclination may be raised to a = w = 69° to pl 
duce overcoverage in the polar regions. The lines~ 
nodes in each case are coincident. 

The principles used for global coverage are readi 
extended to cases of restricted coverage, say, betwe 
an upper latitude \, (small circle B), and a lower la 
tude d; (small circle A). As \, — dz decreases, the 1 
quired number of satellites for both types of patter 
is reduced in general. The photographs (Figs. 9 and 1 
are intended as an aid in the vizualization of types 
pattern coverage. On both, solid white lines indica 
the trace or locus of sub-satellite points on a fix 
Karth. The boundaries of different rings carry differe 
markings. The two patterns shown are directly co: 
parable, i.e., the areas of required coverage are t 
same (A, = 75°, , = 30°). The degree of overlap, A 
may be considered due to the integer constraints on t 
variables of required overcoverage demanded by ott 
factors. 
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Fig. 9. Typical L-type pattern 


Fra. 10. Typical A-type pattern 
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Pattern Stability 


A full treatment of pattern stability is beyond the 
scope of the present paper. One may quite easily 
delineate the cause-effect relationships however. Some 
quantitative conclusions can be made which supplement 
the heuristic value of such an approach. 

As previously defined, pattern stability may be 
related to the time interval between completion of the 
pattern and loss of required coverage due to non- 
Keplerian forces (not accounted for in the model) 
which give rise to differential motion of the ICA’s. 
The definition is extended to include the effect of 
placement inaccuracies. Thus departures from the 
specified or nominal orbit elements during pattern 
attainment may in time cause “holes” to open in the 
coverage area. Table 2 has proved useful in this regard. 


TABLE 2 
ICA Variations 
Effects on ICA 
Cause Positional] 
Size 

Along Orbit Track | Across Orbit Track 
Aa + or —, © a + or —,0 
Ae +or—,P —— + and — 
Aa eee ap ByinGl =. 12 ——— 
AQ == sr Buel ==, 12 ae 
ATa ap OF =, 0 === Sees 
(——— denotes no effect) 


Since the argument of perigee is not defined for a 
circular orbit, we may examine the effect of off-nominal 
values in the other five elements on the ICA position 
and size. These elements are taken to be: 


a, semi-major axis 
e, eccentricity 
a, inclination 
Q, longitude of the ascending node 
To, time of nodal passage from the epoch 


If some element, say a, has a variation Aa, then a 
“1 and —” effect indicates bi-directional changes in 
the ICA position or size irrespective of the sign of Aa. 
A “+ or —” effect depends on the direction of the 
cause. The symbols 0, P, » denote the period of the 
effect. They correspond, in order to a bias shift, an 
orbital periodic variation and a steady drift. One 
notes that only Aa produces a drift effect and that only 
in the ICA position along the track. All others may be 
tolerable due to the intrinsic overcoverage of an 
optimal pattern or “designed out”? by the additional 
overcoverage of an off-optimal pattern. 

In addition to placement inaccuracies, gasdynamic 
drag, higher harmonics in the Earth’s gravitation field 
and third-body perturbations will alter the orbits. It is 
emphasized at this point that orbital changes per se 
do not degrade the coverage capability of a pattern. 
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What is important here is a “differential pertubati 
that acts differently on satellites to change th 
relative positions in the pattern. 

The primary perturbations of a two-body orl 
arising from the Earth’s oblateness are as given 
King-Hele and Merson:? 

(1) secular rotation of the line of nodes 

(2) radial excursions of half orbital period 

(3) dependence of nodal period on the inclination — 

(4) secular rotation of the line apsides 
The last has little significance for nominally circular ¢ 
bits; the third has no effect since ring geometry is inv 
ant with a change in period common to all satellites 
the ring; the second is of quite small magnitude and 2a 
easily be designed out. The regression of the line « 
nodes, 2 = 10 (R/a)*® cos a degs./day, has no effec 
on L-type patterns since all satellites have the same 
and a. The pattern merely “sees” the Earth rota 
more rapidly than once per day if direct orbits an 
employed. The rings of A-type patterns however, wi 
experience a differential nodal shift and in general, 
hole will open in some finite length of time. Som 
preliminary computations indicate lifetimes of th 
order of months for optimal patterns although th 
results are insufficient to serve as a basis for an 
general conclusions. We note this disadvantage ¢ 
\-type as compared with L-type patterns and, at th 
same time, observe that off-optimal A-type pattert 
could increase their lifetime indefinitely. 

A study of the influence of celestial bodies such ¢ 
the Sun and Moon (and perhaps other artificial satellite 
as well) on pattern stability lies in an area of investigs 
tion noted for its difficulty and scarceness of quantits 
tive results. Without attempting to describe the preser 
status of the problem of a multiple-body system unde 
Newtonian forces, one notes that even in the restricte 
three-body case, no complete solution is availabl 
If linearity is assumed so that the separate effects o 
pattern orbits from the presence of other bodies ¢a 
be added, the conditions of the restricted three-bod 
problem are met for each perturbing body in tum 
If, further, the most adverse constellation of the Eartl 
satellite-perturbing body system is supposed, ‘“‘bounds 
may be placed on pattern lifetimes. This approaé 
essentially ignores the dynamics of the problem an 
is useful only in estimating the maximum differenti: 
acceleration two adjacent pattern satellites migl 
experience. 

Gradecak’s work*® would indicate the existence ¢ 
secular perturbations only the node and the argument 
perige—the latter defined on the perturbed orbit wit 
periodic eccentricity variations about the zero poin 
The magnitudes of these drifts are exceedingly smé 
for near-Karth (a/R < 2) satellites amounting to a fe 
minutes of are per year. An additional theoretic resu 
is based on the recognition of “resonant orbits”? whi 
tend to amplify the effect of these small perturbation: 


ris: orbit of semi-major axis a = a, is resonant 
ith, for example, the lunar orbit a = a», when the 
quation va; = yar! has integer solutions (y, v’). 
falues of a, which admit integer solutions are to be 
voided, especially if these solution pairs are composed 
f small integers. These orbits for small » and »’, say 
ach an integer less than ten, are beyond the distances 
ne might expect to attain the limiting patterns of 
ractical significance, namely those of three satellites. 
ower reasonant orbits correspond to very high order 
erms in the perturbation function. These amplified 
ariations would be exceedingly slow. 

The gasdynamic decay of equal-altitude circular 
rbits does not disturb the ICA relative positioning. 
‘part from a consideration of pattern lifetime based 
n the time to re-entry, the tendency to reduce the 
CA size and hence the degree of coverage needs investi- 
ation. Using Peterson’s data’ and values of d ¢/dr for 
= 90° (the worst case) from Fig. 2, the ICA radius 
ecreases on the order of 10—! and 10 degs./yr. for 
Ititudes of 250 and 1000 miles respectively. The 
ffects of these rates are clearly ignorable. 


)perational Considerations 


Congressional testimony® has brought out some of 
he operational complexities of coverage-type satellite 
ystems. Of primary importance is the feasibility of 
lacing a satellite in orbit with prescribed nominal 
alues on all six orbital elements. This, of course, 
elates to the rendezvous (with a hypothetical satellite) 
roblem which has received considerable attention 
ecently. Let us suppose the availability of a single 
vunch site and booster vehicle which give a capability 
f meeting the required semi-major axis, eccentricity 
nd inclination values for a given ring. These elements 
re independent of the time of launch; the other two 
re not. Once each day the site is in position to match 
he nominal nodal longitude of the ring by using the 
ame reference ascent path. The in-ring spacing 
etween two satellites launched by the above procedure 
‘ould vary according to the period of the orbit. Now 
this period divides into the day with a 1/m remainder, 
satellites can be launched on successive days so that 
leir in-ring spacing is uniform. If an attempt fails, 
ne is assured of another opportunity m days later. 

The above method ignores nodal regression. Some 
ispersion in nodal longitude would result if it were 
rictly followed. At satellite altitudes greater then 
ne Earth radius, the nodal dispersion would be less 
1an one degree for each day after the start of the launch 
quence. Although this value does not appear to be 
ary large, a certain amount of overall readiness is 
acessary to establish rings by this method. Nodal 
spersions of the magnitude implied by, say, a week of 
lily launchings can be reduced with quite small 
it-of-plane thrust corrections at the antinodes of the 
bits. These post-placement corrections may be 
wnecessary if the choice of nominal elements produces 


sufficient overcoverage. The number of satellites over 
optimal to produce the type of overcoverage associated 
with nodal dispersions is small compared to that re- 
quired when the in-ring spacing is ignored. Also schemes 
for correcting the latter errors!!! have comparatively 
large velocity and time requirements. It appears, 
therefore, that patterns with the aforementioned 
divisibility property have an operational advantage. 
We will investigate the possibility of meeting this 
condition by considering a numerical example. 

Fig. 2 gives for ¢ = 90°, ¢ = 60° at r = 2.0 Harth 
radi. From Table 1, m is 4 and 5 for optimal L- and 
A-type patterns respectively. The satellite period is 
0.164 days resulting in 0.10 revolution shift per day. 
This shift is less than the 1/m = 0.25 and 0.20 neces- 
sary. An altitude adjustment of less than 100 n.m. will 
produce the proper phasing. This is considered to be 
within an actual design envelope relating altitude to 
ICA. Other examples would, of course, require greater 
percentage adjustments, but the model allows for much 
give-and-take so that a developable preliminary design 
of the system is assured. 

The launch sequencing and ascent guidance problems 
L-type patterns are not as difficult as those for A-type 
since all rings are at the same inclination. The same 
ascent path relative to the ground tracking and guid- 
ance instrumentation may be used. This factor may 
be outweighed by differences in what might be termed 
the intrinsic overcoverage area. As Figs. 9 and 10 show, 
even optimal patterns have regions on the Earth in 
which multiple coverage is maintained continuously. 
This region is in higher latitudes for L-type and 
lower latitudes for \-type. A choice may well depend 
on which region could least easily suffer a loss of cover- 
age in the event of equipment failure. 

The great costs associated with satellite systems of 
the type described leads one to assume that their 
“reason for being” will be of first-rank importance to 
a nation’s political and military position. It is difficult 
to imagine such expenditures tied to other than determi- 
nate and well-defined system requirements and execu- 
tion. One must not take the specification in probabilistic 
terms of, say, an ICBM system accuracy (CEP) as a 
counter-example to the foregoing argument. Here 
perfect execution results in perfect attainment of the 
objective. The existence of CEP-type specifications 
merely implies the realization that perfection is ethereal. 
On the other hand, a probabilistic or so-called random 
pattern degrades the purpose even though the develop- 
ment and execution of the system were thoroughly and 
exactly carried out. 
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Astronavigation—Guidance and 
Control in Space’ 


Robert E. Roberson? 


Abstract 


Within the last few years there has been an increasing 
recognition of the central role of guidance and control 
in astronautics and a marked increase in the number of 
studies of such problems. The bulk of this literature 
sometimes makes it difficult to see how individual efforts 
fit together, but it is generally agreed that there is a 
legitimate discipline—under whatever name—which 
ean be characterized loosely as “space guidance and 
control”. I propose to call this “astronavigation”. The 
purpose of this paper is to dissect out the ingredients of 
astronavigation, to show what problem areas are 
logically encompassed by the whole disciplinary area, 
and how these fit together. The paper is in no way a 
survey of the literature of the field—the examples 
chosen to illustrate certain problems are often taken 
from the relatively early literature on those topics. 
However, the topical breakdown itself provides a con- 
venient framework for such a survey. 


Introduction 


The term ‘astrodynamics” has been applied by 
Herrick [1] to trajectory and orbit analysis of astro- 
nautical systems. Boni [2] coined the term “sidar-me- 
chanics” to apply to the same area, but his terminology 
as not been accepted. This distinction from classi- 
cal celestial mechanics is useful, but it is convenient 
0 recognize six distinct facets: three direct or analysis 
roblems and three inverse or synthesis problems. 


Celestial Mechanics 


The first direct problem is to analyze the motion of a 
ystem under a completely specified set of non-pro- 
ulsive forces following a specified initial state. This 
an be considered to be purely celestial mechanics, 
omewhat enlarged from the traditional discipline to 
clude dissipative and other nongravitational forces. 
It includes selection of basic reference systems (see 
3] and remarks in [1]), examination of the values of 
stronomical and gravitational constants [4], and the 
ethodology of analytical and numerical investigations 
f trajectories and orbits. 

1 Portions of this work were supported by Republic Aviation 


orporation. 
2 Systems Corporation of America—Los Angeles, California. 
3 In truth, the vehicle may be a collection of a number of 
igid bodies and may have many degrees of freedom. However, 
e speak here as if it were a single body. 


Analytical attacks have been made on the effects of 
a number of perturbations which were not of much 
classical interest. Two that have received most (and 
probably disproportionate) attention are the effects of 
planetary oblateness and air drag on the orbit of a close 
satellite. The methods have been largely classical, such 
as variations of parameters and other perturbation 
methods, and the goals normally have been to derive 
the low order secular behavior of the orbit elements. 

Numerical work has been of two general types, 
clearly distinguished by Walters [5] as concerning 
“design (or feasibility) trajectories and precision’”’ 
trajectories. The importance of precision trajectories 
also has been emphasized by Herrick and Baker. The 
present author [6] has commented on the domains of 
special value of numerical and analytical approaches 
and has observed that a fruitful area of further progress 
should be that in which these are combined in a single 
study to exploit the special advantages of each to the 
fullest. 

This area of celestial mechanics also should be under- 
stood to include variations about free-flight trajectories, 
particularly those resulting from changes in initial 
conditions. The literature is replete with developments 
(and redevelopments) of the effects of initial position 
and velocity variations upon the elements of two-body 
orbits or upon certain terminal conditions at a subse- 
quent point on a two-body orbit, topics which have 
direct and important implications to guidance and con- 
trol. There has been less attention to the terminal effects 
of initial-errors following other kinds of trajectories 
(e.g. three-body trajectories), for these must be ob- 
tained by numerical computer studies, but this situa- 
tion has changed rapidly during the past year. Most 
results of this kind are available only as company re- 
ports, but one published example is the earlier work of 
Lieske [7]. 

The celestial mechanics aspects of astrodynamics are 
particularly important in analyzing certain operational 
aspects of space vehicle applications some of which 
bear indirectly upon the problems of guidance and con- 
trol. However, it is in the question of variations about 
or deviations from the nominal that the celestial me- 
chanics type studies impinge most directly on guidance 
and control. 
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The inverse problem that corresponds to the analyti- 
cal area of celestial mechanics is that of using observa- 
tions on orbits or trajectories to determine certain 
parameters in nature, such as Earth’s oblateness and at- 
mospheric properties. This kind of problem has been 
discussed in general terms in the literature but has be- 
come much more pressing with the accumulation of 
data on existing satellites. It seems to be a somewhat 
more difficult problem than originally expected, es- 
pecially when observations are of indifferent quality. 
Extracting “average” properties such as a smoothed 
atmospheric density vs. altitude relationship is possible, 
but details of atmospheric structure (including diurnal 
and latitude variations, atmosphere rotation, and the 
like) and of other dissipative effects are well masked by 
gross air drag and nonconservative perturbations, and 
a forteriori by observational errors. An example of what 
can be done in this problem area when accurate element 
determinations are available is given by recent work of 
King-Hele [8]. There is no need to dwell on this topic, 
but it illustrates a definite region within which astro- 
dynamics and guidance may be relatively unrelated 
problems and corrects the occasional implication that 
astrodynamics 7s space guidance. 


Maneuvers 


The second part of the trajectory and orbit problem 
is concerned with powered maneuvers. Given the basic 
operational concept for a certain astronautical mission, 
what kind of nominal path should be followed, and what 
kind of thrust or steering program is needed so that it 
will be followed under nominal conditions? Here again 
we have the distinction between analysis and synthesis, 
but in this case the analysis part is relatively straight- 
forward and the problem of synthesis tends to be prima- 
rily important. 

The synthesis problem also may be called the prob- 
lem of the ‘‘maneuver concept”. I have elaborated the 
role of this topic in establishing guidance and control 
requirements in [9]. Fundamentally, the formulation of 
a maneuver concept is deciding what kind of maneuver 
sequence is “best”? used to accomplish the mission goals 
(where the system is to go and why) within the limita- 
tions imposed by the system’s physical nature, including 
such factors as rocket engine characteristics and avail- 
able propellant weight. For the case where a high quasi- 
impulsive rocket thrust is used, this means a device of 
impulse directions and magnitudes for a set of thrust 
periods separated by periods of free-flight. 

This situation is more difficult in cases where the 
thrust is not quasi-impulsive, then there are two major 
subcases, one relating to high thrust boost and the 
other to very low thrust space flight where propulsion 
is by ion rocket or similar means. The former is a prob- 
lem which has been treated rather extensively, often 
by methods of the calculus of variations (for example, 
in numerous works of Miele). Although the low-thrust 
case has been approached this way (e.g. in [10]), the 
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choice of a nominal low-thrust steermg program prob 
ably is better done on the basis of resulting guidance anc 
control simplicity than on traditional propulsion o 
time considerations. The author inclines strongly to 
ward this point of view and has heard it well expresseé 
by Levin [11] and others. 

This comment introduces the central difficulty of 
developing a maneuver concept, namely the choice ol 
a criterion of merit. What is meant by “‘best”’ maneuve 
or steering program? Classically, this has been taken 
be a propulsion (characteristic velocity) minimizatio 
Typified by many of the excellent papers of Lawde: 
(e.g. [12]{14]), this kind of criterion is particularly 
appropriate where propulsion capabilities are marginal, 
However, there are many other factors too, such as 
guidance and control practicalities and operational 
considerations, which can become dominant in certain 
cases. For example, even the “simple” case of two-im- 
pulse maneuvers becomes complex when multiple 
“payoff” factors, not necessarily consistent, arise. Ir 
[15] the author discusses the selection of an ascent tra- 
jectory for satellite rendezvous where there is a trade- 
off between characteristic velocity and expected waiting 
time, an operational factor sometimes of extreme im- 
portance. } 

Many other examples can be quoted to illustrate the 
deficiencies of the usual single-variable criteria. How- 
ever, it suffices to say that there is increasing recogni 
tion of the operationsl advantages of certain classically 
‘“non-optimum”’ trajectories, and further developments 
can be expected along these lines. At the moment, the 
case of flight throught the atmosphere probably is best 
under control while that of low-thrust space flight is 
least. 


Low Thrust Maneuvers 


Because of the importance of the low thrust maneu- 
ver problem, it is examined further here with reference 
to existing literature. Consider first the direct problem 
of determining the path when a steering problem is pre- 
scribed. The prescription, in this case, takes the form of 
a rule which relates the orientation of the thrust vector 
to time and to such observables as the instantaneous 
vehicle position and velocity. 

Several basic problems of the direct type have beer 
mentioned time after time in the literature. Ordinarily 
it is assumed that thrust is directed either radially (i.e 
along the line from the vehicle to the gravitating center) 
circumferentially (i.e. perpendicular to above and it 
plane of travel), or tangentially (i.e. along the forwarc 
velocity vector). Assumed steering programs invariably 
are one of these or some combination of these with ¢ 
change from one type to another at chosen times anc 
with thrust levels which are piecewise constant or whiek 
vary in some special fashion which permits closed form 
solutions of the dynamical equations. Past studies have 
had one of two basic goals. Either they have beer 
directed toward finding the path or related characteris 


TABLE 1 
L iterature on the “Direct”? Problem of Low Thrust Trajectories 


Purpose 


Thrust Direction 
q Path, Element Changes and 


| . Related Topics 


Propulsion and 
Related Topics 


Radial 11} 16,17; 18,19) 20; Geet ie 22 
Gircumferential ILE ie Wie 2B ING, i, BRA, Pay Oe 
Tangential 23, 25,-26, 27, 28, 29 29,03) 34107..30 


tics (such as breakaway time, changes in orbit elements) 
or they have concerned some aspect of propulsion (such 
as characteristic velocity). Table 1 summarizes some of 
this literature on the direct problem in terms of these 
characteristics. Some of the references have an ap- 
proach which is purely analytical, others an approach 
which is purely numerical, but most involve a mixture of 
analysis and numerical computation. 

_ Next consider the inverse problem of determining a 
nominal steering program when the origin and destina- 
tion are assigned and when all of the dynamical aspects 
of the motion can be described completely, because, 
there are no unpredictable or error terms in the equa- 
tions of motion. A related problem is to determine the 
steering program which results in specified changes in 
orbit elements. 

There have been three basic approaches to this prob- 
lem. One has been exploratory. A steering program is 
chosen which is “intuitively natural” and a resulting 
trajectory is computed. If the terminal conditions are 
not those desired, the parameters of the program are 
changed slightly. For example, the duration of steering 
in one direction could be prolonged at the expense of a 
subsequent steering mode. Or, the thrust level might be 
readjusted. Then the trajectory is recomputed and the 
terminal conditions are again determined. It is hoped 
in this process that judgment of the investigator will be 
an adequate guide so that the procedure converges upon 
an acceptable trajectory. In the problem of changes in 
orbit elements, the required changes are usually 
assumed small so that analytical approximation meth- 
ods apply. Most of the tentative or preliminary design 
type low-thrust trajectories obtained heretofore seem 
to have resulted from this method. Table 2 tabulates 
the pertinent literature. 

A second method for determining steering program is 
based upon the ingenious investigations of Rodriguez 
(32|. He has observed that at least for the case of trans- 
er to a circular terminal orbit, the initial and terminal 
cinematic states can be represented as single points in 
n energy-momentum phase diagram. It follows that 
y path of states in this phase plane followed by the 
ehicle which contains the initial and terminal point 
is a steering program guaranteed (in the absence of 

eviation-producing errors) to result in a proper tra- 
jectory. Rodriguez shows how steering programs can 
e built up naturally from constant energy and constant 


TABLE 2 
Literature on the “Inverse” Problem of Low Thrust Trajectories 
Bieneuver oF Thrust Direction Sought Re 
Radius change | tangential transfer time vs. | 26 
acceleration 
circumferential | steering mode, 15 
and radial switching point 
Kecentricity intermittent, time to change 26 
change tangential vs. acceleration 
level and 
thrust duration 
radial acceleration level] 18, 19 
Various various, includ- | exploration of 31 
ing some effects 
thrust, magni- 
tude and vari- 
ations 
Mars tangential cir- acceleration vs. | 23 
cumferential time 
Moon tangential reversal points 28 


angular momentum segments of trajectory, together 
with other segments of a “natural” type. (E.g. with 
thrust in the direction of the velocity vector or thrust 
normal to the radial vector between the vehicle and 
gravitating center.) In this procedure it is not possible 
to reduce the problem entirely to graphical and 
analytical means, since certain quantities of interest 
still can be obtained only by machine solution of differ- 
ential equations. However, the Rodriguez method is a 
convenient framework for visualizing the effect of 
changes in the steering program and of systematizing 
the exploratory type investigations. It is felt to be 
worthy of considerable further study, with particular 
emphasis on the representation of the effects of errors 
within the same framework. 

A third procedure might best be called the “rational” 
method, consisting of the generation of a steering pro- 
gram by a purely formal means after setting up the 
problem as one in the calculus of variations. All ap- 
proaches of this class have certain things in common. 
First, one must have an explicit criterion for selecting 
the trajectory. Second, one must establish mathe- 
matically all of the essential constraints on the prob- 
lem. Third, the variational equations, which tend to be 
much too complicated for analytical treatment, must 
be formulated in a way amenable to machine computa- 
tion. Fourth, one must investigate the sensitivity of 
the resulting solution to changes in parameters and see 
whether the maximum or minimum under the criterion 
chosen is a sharp one or a broad one. Fifth, one must 
reconcile the resulting steering program with control 
mechanization practicalities. Published investigations 
of this class seem to be those of Lawdon [13], [14] and 
Faulders [10]. Lawdon uses a minimum energy criterion, 
Faulders a minimum time criterion. The major differ- 
ence otherwise is that Lawdon focuses his attention 
primarily upon the calculus of variations formulation 
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while Faulders gives considerable attention to the 
machine solution problem. 

It is seen from the above commentary on previous 
work that most previous results concern the direct, 
problem of finding trajectories numerically, while some 
investigations have been made of the inverse problem 
determining a steering program given the end con- 
ditions. However, our state of knowledge about the in- 
verse problem is still rather sketchy and considerable 
additional work is required. 


Attitude Dynamics 


Heretofore, ‘‘astrodynamics” has been used ex- 
clusively as referring to the path motion of a space 
vehicle—the motion of its center of mass. However, 
space vehicles actually have six degrees of freedom, the 
additional three degrees of freedom representing rigid- 
body motion about the center of mass.? Dynamical 
problems associated with rotational or attitude motion 
are completely analogous to those discussed above for 
translational or path motion. There is no rational basis 
for failing to include them under the heading ‘‘astro- 
dynamics’’, in spite of the fact that it has not been con- 
ventional to do so. Clearly, trajectory and attitude 
dynamics are intrinsically interrelated during powered 
flight phases, and questions of attitude dynamics are 
implicitly encompassed by the previous area of ‘“‘maneu- 
vers’”’. Table 3 tabulates astrodynamics problems. 


TABLE 3 
Astrodynamics 
Major Areas “Direct”? Problems “Tnverse”’ Problems 
Celestial Me- a) perturbation Inferring trajectory 
chanics analysis by an- characteristics or 
alytical approx- physical parame- 
imation ters from obser- 
(b) numerical com- vation of path 
pilations of 
‘feasibility”’ 
and ‘““precision”’ 
trajectories 
(c) effects of varia- 
tions about free- 
flight path 
Maneuvers (a) determination Development of a 
of trajectory maneuver con- 
given thrust and| cept: 
steering program] (a) impulse se- 
(b) effects of varia- quence 
tions about (b) powered flight 
powered-flight in atmosphere 
path (c) low thrust space 
flight 
Attitude Dynam-| (a) dynamical equa- Inferring parame- 
1¢s tions of motion ters from obser- 
(b) derivation of vation of attitude 
torques 
(c) free-body mo- 
tion analysis 
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On the other hand, the attitude-dynamical analogues 
of celestial mechanics are essentially new to the disci- 
pline of astrodynamics. They could be subsumed under 
“celestial mechanics” by appropriately generalizing 
that term. However, I believe that a new name for this 
area is justified and propose simply “attitude 
dynamics’, understanding by this term free-flight at- 
titude dynamics. 

Problems of attitude dynamics are effectively the 
“foundation problems” of attitude control which have 
been outlined and discussed in [33] and [34]. Such topics 
are included as the choice of suitable attitude reference 
frames, the development of the basic dynamical equa- 
tions of attitude motion (non-trivial if the space vehicle 
contains moving parts), the derivation of torques acting 
on the vehicle from external sources (or internal, other 
than torques deliberately applied for control), and the 
behavior of spin-stabilized satellites. 

These are problems of the “direct”? type, but im 
principle the “inverse” problem exists as well. For ex- 
ample, one might imagine inferring something about 
atmospheric or electromagnetic field parameters or 
other dissipative effects by observing the damping of 
rotational motions. Such damping has been analyzed 
to some extent (see references in [34]), particularly as 
it applies to spinning bodies in the Earth’s magnetic 
field, but principally from the ‘direct’? point of view. 
Probably one reason why the ‘inverse’? problem is 
poorly represented is the observational difficulty 
attached to it, unless the vehicle is specially imstru- 
mented with these effects in mind, compared with ob- 
servations of position. 

In summary, the problems of free-flight attitude 
dynamics represent an important, new, but rapidly 
growing area which properly should be included under 
the heading ‘‘astrodynamics’’. 


Astronavigation 
Terminology 


Even if everything behaves perfectly as predicted, 
there are still guidance and control problems, for the 
vehicle must be made to follow its nominal path with 
its nominal attitude. In fact, though, things do not be- 
have perfectly and deviations from nominal conditions 
are bound to occur. It is not enough to precompute 
steering or torque signals needed to maintain nominal 
conditions—overcoming the effects of deviations is the 
very essence of guidance and control. 

Thus far I have spoken of ‘guidance and control” to 
describe a certain problem area, and the reader prob- 
ably has understood in a general way what was in- 
tended. However, I submit that neither term is quite 
adequate to characterize this field, while both together 
are awkward. What one is really trying to accomplish 
by the operation of “guidance and control” is to cause 
the vehicle to pass through a set of kinematic states 
which satisfy a prescribed set of terminal and perhaps 
en route conditions. A “kinematic state” in free flight is 


described by any set of parameters which permits the 
prediction of subsequent motion, such as initial values 
of position, velocity, orientation and angular velocity. 
Or, iM some cases position and velocity may be replaced 
by a set of orbit elements. We may think of the kine- 
‘matic state as constant during any period of free flight, 
changing during any powered maneuver. (See [13]). 
_ To see that this is more than just “guidance” or 
“control” in the traditional sense, consider Figure 1 
which shows schematically the processes involved in 
carrying the system through a set of kinematic states. 
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Fra. 1. Generalized space vehicle control operations 


This is a generalized diagram in the sense that it repre- 

sents a multi-variable control system. The vehicle is 
acted upon by generalized forces (i.e. forces and torques) 
and its behavior in terms of some set of variables is 
sensed. This may be position and/or velocity relative to 
some convenient coordinate system, orientation or 
orientation rate, or any other physical quantities which 
characterize kinematic states. In one case, these quanti- 
ties are operated upon by generalized filtering and the 
resulting signal is fed to control actuators to stabilize 
the entire closed loop control system. Generally the 
purpose of this during powered flight is to keep the 
vehicle attitude from becoming unstable. This portion 
of the overall control of kinematic state is sometimes 
ealled the “inner loop” or might be called simply the 
“stabilization loop”. All of the portion of the closed 
loop system shown in Figure 1 to the left of the dashed 
line is traditionally the “flight control system”. The 
portion consisting of the actuators, “equalization” 
block, and corresponding sensing operations is some- 
times called the ‘‘flight controller’. Note that for free 
flight attitude control, this loop may constitute the 
entire control system. 

On the other hand, the sensing operation also sup- 
plies outputs which are used for calculating or pre- 
dicting certain “conditions”, such as terminal position 
and/or velocity, orbit elements, or other similar factors. 


These “conditions” are whatever is needed to char- 
acterize the thing that the operation of path control is 
trying to accomplish. Predicted conditions are com- 
pared with nominal conditions and the differences or 
deviations are used as the basis for computing control 
commands. In general the latter consist both of steering 
signals and off-on signals. The control command signals 
are combined with the output of the equalization block 
and the result used as an input to the control actuators. 
The portion of the control system consisting of the pre- 
diction operation, comparison, and generation of path 
control commands is traditionally called the “guidance 
system”’. It is labeled “guidance ‘loop’ ”’ in Figure 1, 
although “loop” strictly is a misnomer—it is not a 
closed loop until it is made so through the control sys- 
tem. 

It is clear from Figure 1 that while there is a rational 
basis for the traditional distinction between guidance 
and control, these really have an inherent unity. They 
are just different loops in an overall control system, and 
the exact point of dichotomy is to some degree 
arbitrary. I believe there are advantages in deempha- 
sizing the traditional division and speaking of the entire 
closed-loop system as simply ‘control’, immediately 
dividing it into “path control” and “attitude control’’. 
This neither abandons the former terminology, nor 
raises problems by extending its meaning in unfamiliar 
ways. The teminology of attitude control is already well 
accepted, and one merely adds “‘path control” as a pre- 
ferred substitute to ‘guidance and (flight) control’. 

Note that it is in the external force, vehicle dynamics 
and nominal conditions blocks, and possibly in the pre- 
diction block, that astrodynamics enters the picture. 
The ‘control’ area, as I have proposed it here, would 
use the astrodynamical results as essential background 
but would be itself concerned with the rationals under 
which the operations in the various blocks are per- 
formed, the physical structure of these blocks, and the 
closed loop characteristics of the entire combination. 

There is the problem of what to call the entire area 
that is neither entirely astrodynamics nor entirely 
control. I propose “astronavigation” as a suggestive 
term which seems suitable despite a possible confusion 
with the interpretation ‘navigation by reference to the 
stars’. It usually begs the question to refer to the dic- 
tionary, but one does find there that “navigate” is to 
direct one’s course or steer, while ‘‘navigation” is the 
act of navigating using calculations of position and 
direction, etc. This is exactly what we do mean by the 
operation we are discussing here whose major ingredi- 
ents are astrodynamics and control. 


Path Control 


The foundations for path control are laid with the 
establishment of a maneuver concept which results in 
a nominal flight path. However, this astrodynamical 
question is only a first step, and beyond it there lie 
five areas of investigation which are unique to the path 
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control problem. The first is a careful examination of 
the nature of the disturbances acting on the vehicle 
to cause it to depart from its programmed trajectory. 
The second is a study of possible methods that might 
be used for sensing its deviation from the desired path 
or its absolute position and velocity relative to a pre- 
scribed reference frame. The third is a consideration of 
the characteristics of the devices that might be used to 
change its path, the actuators. The fourth is the de- 
velopment of a correction philosophy, which is not to 
be confused with the maneuver philosophy chosen in 
the absence of any error considerations. The fifth is a 
performance analysis in terms of the various sources of 
error that may be expected. 

The disturbance analysis is important for several 
reasons. The qualitative nature of the disturbances 
can influence the collective action against them or im- 
pose specifications on the operation of certain parts of 
the system. For example, in a low thrust system one 
may find that major disturbances arise from random 
stochastic variations of line of thrust relative to the 
vehicle frame. The best remedy in this case might be 
to use the disturbing mechanism also as a corrective 
mechanism or to tighten the control of the thrust vector 
to a point where the stochastic effects are negligible. 
The quantitative magnitude of the path deviations re- 
sulting from the disturbances is important both as it 
effects the choice of sensors and a control concept. As an 
example, if no disturbances act to carry the vehicle 
out of its nominal plane of motion it may be possible 
at certain levels of accuracy to dispense with sensors 
for detecting such out-of-plane motion, or at least to 
give them a rudimentary form. If path deviations re- 
sulting from disturbances are small, but the path con- 
trol requirements of the mission are stringent, it is 
necessary to use sensors with a low threshold level, i.e. 
high sensitivity. If deviations are large and accuracy 
requirements are still stringent, it is necessary to use 
sensors with a larger dynamic range and low scale factor 
error. In all, the analysis of disturbances has a number 
of implications to the later choice of sensors. 

In order to illustrate the possible connection between 
deviations and the choice of a correction concept, one 
may anticipate some comments below about correction 
possibilities. An important one is based upon the ex- 
pansion of true position about nominal position in a 
Taylor series in the deviations and basing the correction 
calculations on only the first order terms in small devia- 
tions. Clearly, if the deviations become large the second 
and higher order terms of the expansion become much 
more important and it may introduce serious errors if 
the correction is made only on the basis of the linearized 
form. Large deviations also assume particular im- 
portance for the case of low thrust vehicles, since it 
may become impossible to choose any path control 
concept which permits their correction within the thrust 
limitations of the propulsion system. The transition 
from existence to non-existence of a correction concept 
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is obviously a rigorous one and nicely illustrates that 
the magnitude of deviations indeed can have an 1 
portant effect. Actually, one might expect a whole 
series of gross changes in the required concept as the 
limits of thrust capability are approached. 

Let us turn now to the question of sensors. It is per- 
haps unfortunate that so much stress has been placed 
in the literature on the methods of sensing position and 
velocity of astronautical vehicles. This has resulted in 
tendency to regard this sensing operation alone as “the” 
problem of space navigation, whereas “it is actually 
just one facet whose importance should be neither 
minimized nor over-magnified. Considerable ingenuity 
has gone into the question of establishing a location im 
space. Broadly, proposed methods are of four majot 
classes: sightings on distant bodies (effectively point 
sources), measurements on ambient fields of the vehicle, 
line of sight observations on nearby bodies, and ob- 
servations on the vehicle from other sites. 

Stars may be used as point sources for establishing 
orientation, but position inference by similar observa- 
tions involves distant planets. The process is essentially 
one of triangulation, using one’s knowledge about the 
geometric location of the planets. The literature con- 
tains both the basic equations to be solved and 
expressions for the position errors resulting from angulaz 
measurement errors and the error in the astronomical 
unit of distance. Velocity can be inferred most prace- 
tically by differencing successive position measure- 
ments, but it also has been suggested that the Brewster 
effect or doppler shifts might be used. 

The simplest examples of position and velocity de 
termination by field sensing occur in the neighborhood 
of planets where an atmosphere and magnetic field 
might exist. Even the gravitational field is of potential 
use to infer distance from a nearby body if sufficiently 
precise measurement of the gravity gradient can be 
made aboard the vehicle, e.g. by pairs of very sensitive 
accelerometers. Still another method is to determine 
distance from a radiating body (typically the sun) by 
absolute measurement of received radiation at the ve- 
hicle. 

Line of sight observations on a close body offer a 
number of possibilities, such as triangulation by either 
angle or absolute distance (by radio) measurements 
using three points or sources whose geographic loca- 
tions on the body are known. Observations on the pe- 
riphery of the presented disc, in the manner of horizor 
scanning often discussed for attitude determination 
offer the possibility of determining radial distance pro: 
vided the diameter of the body is known and one car 
accept such systematic errors as those from oblateness 
and refraction. 

The final class, observations on the vehicle from othe 
locations, is a particularly important one, A large num: 
ber of possible methods is represented, including optica 
and radio methods. More than one external locatior 
may be involved, perhaps the planet of origin or des 


, 
- or on another space vehicle. Cooperative ar- 
ngements using beacons or transponders on the ve- 
hicle would be included here. It should be recognized 
that if these sensing methods are to be included in the 
control loop, an additional communication of command 
information to the vehicle is required. This, among 
other operational considerations, would suggest that 
such sensing is of value principally in low frequency 
applications such as occasional impulsive course cor- 
rection. But frequent position and velocity information 
generally is necessary only to initiate maneuvers near a 
planet, so periodic position checks supplemented by 
interim prediction may suffice in any case. One 
advantage of this class is that complete geographic 
location including not only radial distance but also the 
equivalent of a latitude and longitude can be obtained, 
information which cannot be extracted from certain of 
the other methods. 
One other point about sensing also must be remem- 
bered, namely that in path control the sensing plays 
two roles. It not only determines position and/or veloc- 
ity, but also establishes the reference frame with respect 
to which the path control operations should be con- 
ducted. For the latter, there is a possibility of inte- 
grating path and attitude control sensors as shown in 
Figure 1. 

- The sensors discussed above form the portion of the 
path control problem traditionally called guidance. In 
the same way, the actuators together with some of the 
control loops comprise the traditional flight control 
area. Actuator characteristics are of obvious importance 
in making any error studies of a complete system, but 
also may have an effect upon the choice of a control 
concept. For example, the relative accuracy with which 
a thrust device can be pointed and turned off or on 
could determine which of these two operations was 
stressed in this concept. There are two basic forms 
which actuation takes in path control: one is to change 
the direction of the line of thrust so as to change the 
direction in which the center of mass of the system 
moves, and the other is to initiate and terminate thrust. 
The direction of the thrust in space normally is changed 
by changing the orientation of the entire vehicle frame: 
this, in turn, is accomplished by changing the direction 
of the line of thrust relative to vehicle fixed axes so that 
it no longer passes through the vehicle center of mass. 

he physical mechanism for accomplishing the change 

f direction relative to the vehicle can be either a gim- 

alled engine or fixed side thrust jets. The fact that the 

otational mode of the vehicle arises in this problem im- 

lies that there are really two control loops to consider. 

ne of these is the relatively high frequency stabiliza- 
ion of the vehicle attitude during powered flight, so 
hat. the line of thrust is kept properly oriented in 
nertial space. The other is the lower frequency guidance 
ode which pertains to the gross motion of the center 

f mass. Both of these are proper ingredients of the 

ath control problem. 


The path control concepts deserve considerably more 
attention than can be given them in this brief descrip- 
tive review. It is convenient to visualize three varieties: 
concepts for the establishment of a prescribed initial 
state (e.g. position and velocity) where the path in- 
cludes a main boost period, those for the change of 
state by means of impulsive or quasi-impulsive maneu- 
vers; those applicable to low thrust correction systems. 
In the first, the concept for the correction of deviations 
from a nominal path may be any of those applicable 
to a ballistic missile, which is the case that probably 
has received the most extensive past study. In the open 
literature, the concept generally described is the so- 
called “delta method”. The position and velocity devia- 
tions from the nominal trajectory are monitored during 
main boost. If cut-off should occur at any instant, the 
corresponding deviators at the terminal point (i.e. a 
point following a coasting period at which the mission 
ends or some further maneuver is initiated) are given a 
Taylor series in the instantaneous position and velocity 
deviations, often approximated closely by just the 
linear terms in the series. Predicted deviations at the 
terminal point can be used in at least two ways. For 
instance, the steering can be modified (e.g. by changing 
the line of thrust as some linear function of observed 
deviations) and the cut-off time can be advanced or 
retarded so as to make at least some of the terminal 
deviations equal to zero. Alternately, the original 
steering and cut-off program are unchanged while the 
observed cut-off deviations are used to compute a new 
course of action for the maneuver to be performed at the 
terminal point. In general, the actual concepts for a 
specific mission involve a combination of these. 

At the same time, high frequency stability of rota- 
tional modes is necessary during boost, so it 1s necessary 
to include an additional control loop which assures that 
the vehicle attitude is stable relative to the value pre- 
scribed by the course of action above. It is well known 
that this is an exceedingly difficult task because of ve- 
hicle flexibility and change in its parameters as fuel is 
used, but obviously it is a problem which has been 
solved in at least some cases. 

The change of elements by impulsive maneuvers is an 
elementary topic in astrodynamics and it would appear 
that the path control concept is quite straightforward. 
This is largely true: one need merely choose an initial 
vector velocity increment which makes the desired 
changes. The first difficulty is that only three velocity 
components are subject to control, whereas it may be 
desired to change more than three elements. The second 
is that any velocity component generally changes ele- 
ments other than those desired. These considerations 
mean that more than one velocity impulse generally is 
needed, and that several may be desirable. Errors in the 
physical realization of any impulse affect the state of 
the system at the time the next impulse normally would 
be delivered. Thus the next impulse plays not only its 
nominal role (a purely astrodynamical question) but 
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must also correct the errors resulting from previous 
maneuvers. This, then, is the problem of the maneuver 
concept from the path control point of view. 

Very little has been said in the literature about this 
propagation of errors through a general maneuver 
chain, and even less about correction philosophies. The 
case of low thrust error correction also is a topic which 
has received only a little attention. 


Attitude Control 


I have dwelt on the attitude control problem area in 
considerable detail elsewhere [34], so that only a brief 
review is appropriate here. I have found it convenient 
to divide the topic into the following areas: funda- 
mentals of foundation problems, disturbing torques, 
attitude sensors, actuators, and control system con- 
cepts for controller synthesis. In many ways these 
closely parallel the major areas into which path control 
has been divided above. 

The foundation questions involve primarily the con- 
trolled element dynamics, including the choice of an 
operationally suitable reference system. They also en- 
compass performance measures and synthesis criteria. 
The development of the attitude dynamics in the form 
of differential equations of attitude motion is not a 
trivial matter since the vehicle ordinarily consists of a 
collection of parts that can move relative to one 
another. It should be emphasized, though, that even 
if the vehicle were a single rigid body the controlled 
element dynamics could not be idealized simply as a 
transfer function—the full differential equations gen- 
erally must be used. Two major varieties of attitude 
reference are common: the “body-centered” system in 
which one axis is directed approximately toward the 
center of a nearby body and another axis is in the for- 
ward direction of motion about this body; and the in- 
ertial system in which the orientation of the reference 
frame is fixed in inertial space. Within the former class 
a variety of detailed possibilities exist. The remaining 
foundation problems of performance representation 
and synthesis criteria present problems which have not 
yet been treated adequately in connection with attitude 
control systems. 

Disturbing torques on satellite or space vehicle arise 
from many physical effects. For a vehicle containing 
moving parts, especially crew members, this can be ex- 
pected to be the major source of disturbance. Other 
torques arise from non-uniform rotation of the orienta- 
tion reference frame, gravity gradient effects, electro- 
magnetic interaction with ambient fields, aerodynamic 
drag in the fringes of an atmosphere, bombardment by 
particles or radiation (especially solar radiation). The 
analysis of the disturbing torques has been discussed 
in the previous reference and it is wished to emphasize 
only one point here. This is the fact that the torque 
structure on the vehicle is extremely sensitive to the 
configuration design and to such operational details 
as the regions of space through which the vehicle passes. 
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One would expect the dependence on mission but migh 
not be aware of the importance of vehicle design detail 
Considerable attention has been given to the analysi 
and development of sensors suitable for determinin 
attitude, particularly attitude relative to a close bod. 
when the attitude reference is to be body-centerec 
However, one need not consider all of the resultin 
methods since one of them, the horizon scanner, seem 
to suffice to determine the vertical for most body 
centered applications. Even for this one type of methoc 
various mechanizations are possible, but. the propert; 
they have in common is that they somehow determin 
the center of the apparent disc of the nearby body as i 
is presented to an observer or instrument in the vehicle 
Complete attitude determination includes yaw as we 
as vertical, so the horizon scanner or equivalent devie 
must be supplemented by a yaw sensor, one type 0 
which is a gyrocompass. Rate sensing sometimes is im 
portant and can be done by rate gyros or celestial rat 
sensors, the latter, based on the apparent drift of th 
star field, appearing particularly attractive from thei 
overall design characteristics. Alternatively, one ca 
dispense with rate sensing and attain the necessar 
control system damping by differentiation of attitud 
signals within the controller. If the desired attitud 
reference is inertial rather than body centered, th 
most attractive sensors appear to be star trackers 
possibly supplemented by an inertial platform. 

In general, three kinds of control torques need to b 
applied. First, one must apply relatively high torqu 
levels to stop any initial tumbling of the vehicle. Seconé 
one must apply torques of somewhat less magnitud 
and generally of oscillatory character to counteract th 
disturbances during the operational life of the vehicle 
These two cases would suffice if there were no constan 
torque components acting over a long period of time 
but since there normally are such torques the tota 
angular momentum storage capacity of any control de 
vice would eventually become saturated. This means 
in effect, that a third type of torque need be applied— 
“trimming” torque to prevent angular momentum 
saturation of the primary actuator by providing a pat 
for excess angular momentum to escape from the vehicl 
to inertial space. 

The problem of controller synthesis is. the leas 
studied of the major problems of attitude control, pri 
marily because of the complexity of the controlle 
element dynamics and the fact that it is not yet know 
with what degree of approximation one can represen 
the vehicle by a conventional transfer function. How 
ever, some synthesis techniques are in the process ¢ 
development. Analog computer studies of control sys 
tem response using various controller functions hay 
been made by a number of organizations. From exten 
results it appears that controller design is not @ 
formidable a problem as the present author had orig 
nally envisioned, although it is not a trivial one. On 
can design the system so that its accuracy is essential 


iat, of the sensors provided enough control power can 
ye made available. However, in practical systems, it 
en is preferable to relax the accuracy requirements in 
wor of the considerable saving in power consumption 
that results from the broader accuracy tolerances. Thus 
che most economical system from the point of view of 
power consumption (and space vehicles do tend to be 
power starved), is one which permits a slow limit cycle 
oscillation through several degrees. On the other hand, 
uecuracies of a fraction of a second of arc seem to be 
attainable if one is willing to pay the price in power. 
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ANNOUNCEMENT 


The Journal of the Astronautical Sciences wel- 
comes papers on any aspect of astronautics. Con- 
tributions should be original, generally quantita- 
tive in nature, and should satisfy high standards of 
scholarly excellence. Papers are especially solicited 
in the fields of space flight mechanics, space vehicle 
design, space physics, propulsion, guidance and 
control, communication, space medicine and astro- 
biology, and applications of astronautical systems. 
However, any other papers concerned with astro- 
nautical investigations will be considered. 


The Journal of the Astronautical Sciences 95 


Technical Notes 


Optimum Two Impulse Ascents 
to Circular Orbits 


J. P. Carstens* and T. N. Edelbaum* 


Abstract 


The problem of launching an earth satellite into a 
circular orbit of arbitrary radius and angle of inclina- 
tion with respect to the launch point has recently been 
posed by Wolfe and DeBra [1]. It is the intention of this 
note to indicate that, contrary to the statements in [1], 
target orbit interception at the apogee of the transfer 
ellipse is not generally optimum, and leads to non- 
minimum values of required total characteristic velocity. 


Discussion 


Releasing the restriction that interception of the target or- 
bit must occur at the apogee of the transfer ellipse yields the 
results shown in Fig. 1, which is compared to similar data 
presented in Fig. 5 of [1]. Figure 1 was determined numerically 
on a high speed digital computer by allowing both @ and B, 
to vary for any given combination of 7 and ¢min. The assump- 
tion of orbit injection at the apogee of the transfer ellipse is 
evidently justified only for very low values of tmin (exactly 
so for imin = 0), and for all values of ¢min with large r. 

The values of @ for the optimum maneuver are compared 
to the values required for injection into orbit at apogee in 
Fig. 2. @opt drops below 90 deg for small values of r even 
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for apogee injection, in variance with Fig. 4 of [1]. Optimur 
values of a second trajectory variable, 62, the angle betwee: 
the velocity vector in the transfer ellipse and the radiu 
vector at the point of interception with the target orbit, ar 
presented in Fig. 3. B2 < 90 deg corresponds to the satelliti 
vehicle being lofted to a greater altitude than the targe 
orbit, the second impulse being given at the descendin 
node. For injection at apogee, B» = 90 deg. 

The maximum difference between the two ascents evidentl 
occurs at r = 1.0 and tmin = 90 deg. For this case a savin: 
in normalized characteristic velocity of 0.152 is availab 
with the optimum maneuver over the trajectory of [1]. Fe 
low altitude earth orbits this would correspond to a cha 
acteristic velocity saving of about 4000 fps and would allo 
typical increases in payload with chemical rockets of fror 
40 to 55% [2]. 

A more complete discussion of this problem and the opt: 
mum transfer ellipses may be found in [3]. : 


Nomenclature 


tmin Minimum possible angle between plan containin 
launch point and center of earth and plane of targe 
orbit. 

i Nondimensional radius, normalized with respect t 


radius of launch point. 

Vo Cireular velocity at launch point 

AV Total change in velocity of first and second impulse 

B, Angle between velocity vector and radius vector ¢ 
launch point 

B., Angle between velocity vector and radius vector 
interception of target orbit 

6 Angle between radius vector of launch point and radi 
vector at point of interception of the target orbit. 


ie a a 
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Comment by D. DeBra* on the 

Note by Carstens and Edelbaum: 

Optimum Two Impulse Ascents to 
Circular Orbits 


When Mr. Edelbaum sent us a copy of this note he pointed 
out that the case where imin is equal to 90 degrees and r is 
equal to 1.0 can be easily verified as a counter example to 
Fig. 4 and Fig. 5 of [1]. In this case 7 must be constant and 
the minimum characteristic velocity occurs for the minimum 
semi-major axis (energy) of the transfer orbit. The maximum 
error in our results also occurs for this case as pointed out 
in this note by Carstens and Edelbaum. 

Our analysis of the minimum characteristic transfer orbit 
was checked for a value of r larger than 3.0. It can be seen 
from Fig. 1 why this check did not help to discover the error 
since the two results rapidly converge above this value of the 
radius. All of which goes to show that ‘the usual methods 
for finding extrema” must always be used with caution and 
checked at many points to avoid missing errors of this type. 
Our analysis was performed with the additional constraint 
of second maneuver at apoapsis and therefore led to non- 
optimum transfer orbit as reported here. 

It is surprising that as r became large the two results should 
apparently coincide so that the asymptotes of Fig. 4 of [1] 
are still correct. 

The results of Carstens and Edelbaum are interesting, and 
for small 7, an important correction of [1]. 
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Nov. 21-26 Colloquium on Space Research—Argentine National Committee for Space 
Research and Argentine Interplanetary Society, Buenos Aires, Argentina. 


Dec. 5-8 
ington, D. C. 
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1961 
Jan. 16-18 
Texas. 
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ington, D. C. 


ARS Annual Meeting and Astronautical Exposition, Shoreham Hotel, Wash- 


Symposium on Lunar Exploration—Sponsored by American Astronautical 
Society. To be held as part of the 127th Annual Meeting of the American 
Association for the Advancement of Science, New York, N.Y. 


7th Annual Meeting, American Astronautical Society, Sheraton Hotel, Dallas, 


AAS Special Symposium on Space Vehicles in lonized Atmosphere, Wash- 
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NASA TO CO-SPONSOR 
LUNAR SYMPOSIUM 


The Lunar Exploration Symposium 
to be held by the AAS as part of the 
127th Annual Meeting of the American 
Association for the Advancement of 
Science will be co-sponsored by NASA. 
Ross Fleisig, Chairman of the Sympo- 
sium, has announced the meetings will be 
held on December 27, 1960 in the Bowman 
Room of the Biltmore Hotel in New York 
City. In addition to morning and after- 
noon sessions, an evening panel discus- 
sion on Bio-Astronautics will be held. 
Chairman of the Panel will be Alfred M. 
Mayo, Assistant Director for Bio-Engi- 
neering in the NASA Office of Life 
Sciences. 

The preliminary program for the Sym- 
posium is as follows: 


December 27, 1960 


SESSION | - LUNAR 
EXPLORATION, 9 A.M. 


CHAIRMAN: Dr. Hugh L. Dryden, Dep- 
uty Administrator, National Aeronau- 
tics and Space Administration 


PAPER | 


Scientific Objective of Lunar Explora- 
tion, ROBERT JASTROW, Goddard 
Space Flight Center, NASA 


PAPER II 


Distribution of Interplanetary Dust in 
Cislunar Space, 8S. F. SINGER, De- 
partment of Physics, University of 
Maryland 


PAPER Ill 


Comparison of Special Pertubation 
Methods in Celestial Mechanics with 
Special Application to Lunar Orbits, 
SAM PINES, Applied Mathematics 
Section, Republic Aviation Corpora- 
tion 


CHAS Uattoual 


SOCIETY 
NEWS 


PAPER IV 


Radiation Shielding for Lunar Space- 
craft, T. G. BARNES, A. L. BARA- 
ZOTTI and E. M. FINKELMAN, 
Grumman Aircraft Engineering Cor- 
poration 


SESSION II - LUNAR SPACE 
SPACE SYSTEMS, 2 P.M. 


CHAIRMAN: Robert Young, Executive 
Vice President; ACF Electronics Di- 
vision, ACF Industries, Inc. 


PAPER | 


A Family of Radio Isotope-Fueled 
Auxiliary Power Supplies for Lunar 
Exploration, JUSTIN L. BLOOM, 
Martin Company 


PAPER Il 


Extending the Range of Radar Beacon 
Tracking for Lunar Probes, N. 8. 
GREENBERG, ACF Electronies Di- 
vision, ACF Industries, Ine. 


PAPER Ill 


Horizon Trackers for Lunar Guidance 
and Control Systems, K. H. KUHN, 
Sperry Gyroscope 


PAPER IV 


Lunar Soft Landing Guidance Sensors, 
GORDON BURTON and ARTHUR 
BARABUSH, Raytheon Missile Sys- 
tem Division, Raytheon 


PANEL DISCUSSION 
Is There a Need for a Manned Space 
Laboratory?—December 27, 1960, 8 P.M. 
MODERATOR—A. M. Mayo, Assistant 
Director for Bio-Engineering, NASA 
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SEVENTH ANNUAL 
MEETING TAKING SHAPE 


& 


PAPERS GIVE PROMISE OF 
TOP-FLIGHT MEETING : 


The call for papers for the Sevent 
Annual Meeting to be held at Dalla 
January 16-18 went out on July 7, an 
the excellent response indicates AA 
can expect sessions of unusual interest 


LEADING AUTHORITIES SELECTED 
AS SESSION CHAIRMEN 


The Southwest Section, host for th 
meeting, is striving to make the Sevent 
Annual Meeting a session of exceptione 
benefit to the aerospace industry. As a 
indication of this, a number of the ng 
tion’s leading authorities in the field ¢ 
astronauties have accepted posts as sé: 
sion chairmen to assist in obtaining to 
quality papers covering virtually ever 
phase of space technology. The AAS | 
proud to announce these men as chairme 
of the technical sessions: 


Dr. Samuel Herrick, University ¢ 
California, Astrodynamics Session 


Dr. Otto Struve, National Radio As 
tronomy Observatory, Cosmolog 
Session 


Dr. 8. F. Singer, University of Mary 
land, Solar System Physics Session 


Dr. Homer E. Newell, National Ae 
onautics and Space Administratior 
Geophysical Exploration Session 


Dr. Thomas Gold, Cornell University 
Lunar and Planetary Exploratio 
Session 


Mr. George H. Stoner, Boeing An 
plane Company, Space System 
Appleations Session 


Mr. Henri Busignies, Internations 
Telephone and Telegraph Corporse 
tion, Space Communications Sessio 

Brig. Gen. Don Flickinger, Air Re 


search and Development Commanc 
Bio-Astronauties Session 


Mr. H. Julian Allen, NASA Ames Re 
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Greek Alphabet 


A a_ alpha (a) N vy nu (n) 
B B beta (b) oe ea (x) 
T y gamma (g) O o  omicron (6) 
A 6. delta (d) I wz pi (p) 
E e¢« epsilon (8) P p_ rho (r) 
Z ¢ zeta (z) = os sigma  (s) 
H neta (8) T 7 tau (t) 
6 6 theta (th) T upsilon (u) 
I ce iota (i) © ¢y phi (ph) 
K x kappa (k) X x chi (ch) 
A 2d. lambda (1) W y_opsi (ps) 
Mwy mu (m) Q w omega (5) 


e (epsilon) = strain 

g,s (sigma) = stress 

7 (tau) = shear stress 

» (mu) = micro 

uw (mu) = Poisson’s Ratio 
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Fractions in the body of the text and fractions occur- 
ng in the numerators or denominators of fractions 
should be written with the solidus. Thus: 


cos (1x/2b) 
cos (ra/2b) 


is the preferred usage. 

The intended grouping of handwritten formulas can 
be made clear by slight variations in spacing, but this 
procedure is not acceptable in printed formulas. To 
avoid misunderstanding, the order of symbols should 
therefore be carefully considered. Thus 


(a + bx) cos t is preferable to cos t (a + bx) 


In handwritten formulas the size of parentheses, 
brackets and braces can vary more widely than in 
print. Particular attention should therefore be paid to 
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{la + (b + ex)"] cos ky}? 
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